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REMARKS ON TEE STABILITY OF PLANETARY 

SYSTEMS. 



BY G. W. HILL, ESQ., NYAOK TURNPIKE, S. Y. 

As, in some quarters, quite erroneous views seem to be entertained 
regarding the conditions necessary for tiie stability of the solar system, it 
may be of service to note here, in brief, what is known on this subject. 

It is remai'kable that, although the meaning of stability in statics is 
well known, no one, so far as I know, has ever given a rigorous definition 
of this term as used in dynamics. As applied to the solar system, the 
sense-attributed to it in general seems to involve the idea that the mean 
distances, eccentricities and mutual inclinations of the planets should 
always be comprised within narrow limits. But if this be the proper 
meaning of the word, one is tempted to ask — how narrow? It is plain, 
when we consider the matter more closely, that the distinction between 
stability and instability is one of kind and not of degree. There mast be 
a sharp line separating stable systems from unstable. 

In the first place we must discriminate between two possible significa- 
tions of the term; a system may be stable or unstable with reference to 
the action of foreign forces, or with reference to the mutual action of its 
parts. A slight disturbance from without may cause in a moving system 
only trifling deviations from the previous paths of motion, or the efifect 
may be a greater and still greater departure from them. This is quite 
analogous to the stability and instability of statics. But the stability of a 
planetary system, with reference to its own action, must be defined in a 
way quite peculiar. 

A planetary system is stable when finite superior and inferior limits can 
be assigned to all the distances of the bodies composing it, and that, no 
matter bow long the motion may be prolonged; but if to some or all the 
distances, no superior limit other than infinity, or no inferior limit other 
than zero, can be assigned, the system is unstable. 
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Hence in a stable system, there can be no collision and no indefinite 
separation of the bodies composing it. With this definition of stability 
we see that, in the problem of two bodies, motion in an ellipse is stable, 
but motion in a riglit line, or parabola, or hyperbola is unstable. 

With regard to stable systems, we may enunciate the following prop- 
osition: The coordinates of all the bodies in a stable system, or any 
functions of them which remain always finite and continuous, can be devel- 
oped in infinite converging series of periodic terms, each of the form 
KM,(k z! + /?;; lu^ and /? being absoktely constant; and the argument 
h t -\- ^'m always composed, as a linear function with positive or negative 
integral coefficients, of other arguments, whose number never exceeds 
3w. — 3, w being the number of bodies iu the system. 

With regard to the convergence of these series, it must be understood 
that it is asserted only when K and Ic are taken as wholes ; if these quan- 
tities are expressed in infinite series involving the powers and products of 
certain parameters, these series may cease to be convergent long before 
the system passes from stability to instability. But the convergence of 
these is altogether another question. 

If mathematicians had succeeded in completely integrating in finite 
terms the diiferential equations of motion of a system of material points 
acting on each other in accordance with the law of gravitation, the con- 
ditions under which the motion is stable or unstable could be immediately 
assigned. But there is scarcely any reason to expect that this will ever 
be accomplished, and perhaps it is in the power of analysis, at present 
even, to demonstrate its imposibility. For, just as, from the fact that the 
equation, sin a; = 0, has an infinite number of roots, it may be confidently 
asserted that sin x cannot be represented in finite terms by algebraic 
functions; or that, because an elliptic function possesses the property of 
double periodicity, it cannot be equivalent to any finite expression involv- 
ing circular or logarithmic functions; so it is probable that the functions, 
defined by the differential equations of planetary motion, have properties 
that cannot belong to any finite expression involving quadratures. If this 
should be the case, all attempts to arrive at a complete solution, in finite 
terms, of the famous problem of three bodies, must prove as abortive as 
those made to square the circle or to express elliptic integrals in circular 
and logarithmic functions. 

The solution of the general problem, given the initial positions and 
velocities of a system of material points, to determine whether the ensu- 
ing motion is stable or unstable, in the sense we have attribi;ted to these 
words, does not seem to have engaged the attention of geometers. It does 
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not however demand tlie complete integration of the differential equations 
of motion. Thus where the system contains two material points only two 
integrals are needed, that of the conservation of living forces, and that of 
the conservation of areas. And the criterion of stability is here exceed- 
ingly simple; it depends on the values of the constants annexed to com- 
plete these integrals; to insure stability the first must be negative and the 
second must not vanish. 

The question of stability is intimately connected with the values of the 
coefficients k in the periodic terms of the series into which the coordinates 
can be developed in the case of stable motion. If we substitute for the 
coordinates in the differential equations indeterminate series of this form, 
we shall arrive at a number of equations exactly sufficient to determine 
all the coefficients iT and the 3 « — 3 independent ^'s in terms of Sw — 3 
arbitrary constants which must be given by the initial state of the system. 
Eliminating all the coefficients K, we shall have left 3 w — 3 equations 
determinining the 3 w — 3 independent ^'s. It can be readily shown that 
in these equations Ic appears only in even powers. Then if the initial 
positions and velocities are such that they make these equations afford 
positive and finite values for all the quantities W, and if it be granted that 
these are the proper roots to take, motion of the system is undoubtedly 
stable. But if these equations afford negative values for some or all of 
these quantities, and it be granted that these are the proper roots to take, 
the system is necessarily unstable. In the case where some of the quan- 
tities W vanish, the question of stability must be determined from other 
considerations. Whether in any case it is proper to take the imaginary 
roots of these equations for the quantities A^, or whether, like some anal- 
ogous equations in the theory of heat, they have no roots of this kind, is 
a point which is not yet clear. In all this, as the equations determining 
the ^'s can only be obtained in the form of infinite series, it must be shown 
beforehand that the constants, defining the initial positions and velocities 
of the bodies of the system, enter into these equations in such a way that 
they render the series convergent; else any conclusions, as to the values 
of the /?;'s, deduced from them, are not legitimately established. 

The number of conditions, necessary to insure stability, of course 
increases with the number of bodies composing the system. In all cases 
the constant annexed to the integral of living forces in relative motion 
must be negative. It is needless to say that in the solar system this con- 
dition is fulfilled. Certain popular writers have got it that incommensura- 
bility of mean motions is a sine qua, non of stability; but I am not aware 
that this has been asserted by any geometer or astronomer of note. 
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This mistake doubtless arose from noticing that the near approach of 
mean motions to commensurability produces inequalities having large 
coefficients through the division by the small mean motion of the argu- 
ment. But it will not do to assume that these coefficients increase beyond 
every limit when the mean motion of the argument diminishes without 
limit, or that when this vanishes, there are terms in the planetary elements 
proportional to the time. 

Let us illustrate this point more at length. If a is the mean distance 
of one of the planets, and d an argument whose mean motion nearly or 
altogether vanishes ; then so far as this argument is concerned we may 
have the equation 

da A • a 

-—■ — A sin d. 
d t 

If the mean motion of Q is supposed to vanish, the integral of this equa- 
tion is often written, a = a^ -\- A %xa. d . t^ a^ being the value of a at the 
origin of time. But although this treatment is allowable when we wish 
to find the value of a for small values of t, it will not answer when the 
object is to discover whether a is a periodic function of t or not. For it 
has been assumed that A is constant, whereas it is a function, not only of 
a, but of all the other varying elements which define the dimensions of 
the orbits ; also d is not constant, although its mean motion vanishes, for 
its periodic inequalities may have some eifect here. If the approximation 
were carried further, it would be found that there were terms in a multi- 
plied by ^^, f, &c., and that thus it would be more exact to write 

a = a^ -^ Ami d . t Ar B f -^ Gt^ +.. 

What if the right member of this equation should turn out to be the 
development of a periodic function of ^? This in fact is the result in a 
large class of cases. Thus it is plain that the equation, 

da A . n 
-i— = A sm c/, 
dt ' 

must be treated as a differential equation; that is its right member must 
be regarded as an unknown function of t as well as its left. 

But several instances in the solar system of commensurability of mean 
motions, without resultant instability, ought to have prevented this mis- 
take. The three inner satellites of Jupiter have mean motions generally 
granted to be exactly commensurable, yet the system is not supposed to 
be unstable. Again what prevents our regarding the moon as a planet 
revolving around the sun, and our attributing its being sometimes in ad- 
vance sometimes behind the earth, its having a radius vector, sometimes 
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greater sometimes less than that of the earth to the perturbing influ- 
ence of the latter? In this view the period of revolution of the moon 
about the sun is precisely equal to that of the earth; yet there is no insta- 
bility here. If the planet Venus were moved outward from the sun, until 
its mean distance from this body became nearly equal to that of the earth, 
and if at the same time their eccentricities and longitudes of perihelia 
were so nearly equal as to permit their being for some time in the vicinity 
of each other, the effect of their mutual action would be to make the mean 
values of these elements rigorously equal, and each planet would become 
a satellite to the other. Instability would not result from this disposition. 
There are moreover two remarkable particular solutions of the problem 
of three bodies, in both of which the periods of revolution of the two 
planets are exactly equal, without instability ensuing. These solutions 
have been developed by Laplace (M^canique Celeste, Book X, Chap. VI)- 
The first, in its stable form, may be stated thus : — Two planets may move 
in the same direction about the sun in two equal ellipses lying in the same 
plane, having their foci at the center of the sun, and their greater axes in- 
clined at an angle of 60°, provided the planets are at the same time in 
corresponding points of their orbits, so that they together with the sun are 
always at the vertices of an equilateral triangle. The laws of motion are 
the same as in the case of two bodies, but the common mean motion is 

given by the formula, 

^ __ \m -\- m' -{- m" 
\ a? 

where m, m', m," are the three masses, and a the common mean distance. 

The second is stated thus: — Two planets may revolve in the same di- 
rection about the sun, in similar ellipses having their foci at the center of 
this body, and their greater axes coincident in direction, provided that the 
ratio of the axes is determined by a root of a certain equation of the fifth 
degree, involving in its coefficients the ratios of the three masses. The 
planets must be at corresponding points in their orbits at the same time, 
so that they and the sun always lie in a right line. The laws of motion 
are the same as for two bodies, but the common mean motion is given by 
a complex expression involving the root of the equation just mentioned. 

It may be noticed that Lionville has shown that the latter of these solu- 
tions is unstable in the sense we first attributed to this word, that is with 
reference to slight disturbances from without.* It is probable tliat the 
first solution is in the same case, but I do not know that it has ever been 
discussed in this respect. 

*See Conn. <Xe$ Temps for 1845, or LionvlUes Journul, First Series, Vol. VII, p. 110. 
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Let us now establish in a clearer light the fact that the commensura- 
bility of mean motions does not necessarily produce instability. The so- 
lution of the problem of three bodies can be reduced to the integration of 
a system of eight differential equations of the first order; and by a suit- 
able selection of variables, these may be made to take the canonical 
shape; that is the differential coefficient of each varying element, with 
respect to the time, will be equal to the positive or the negative of the par- 
tial differential coefficient, with respect to the conjugate element, of a 
function B, analogous to, but not identical with, the disturbing function in 
perturbations. Thus the eight elements are divided into two classes; 
four being functions of the mean distances and eccentricities, relate to the 
dimensions of the two orbits; while the other four, their conjugates, are 
simply the elementary arguments of the periodic terms contained by B in 
its developed form. The selection of these last may be made arbitrarily. 
If wc take one of them, as 6, to coincide with an argument of B^ whose 
mean motion nearly or exactly vanishes, and call the element, conjugate 
to this, 0, we shall have the two differential equations 

d_e^dB dd ^ _ dB 

dt dd' dt d &' 

Let us now suppose that B is reduced to its terms which have only 6 
in their arguments ; then 

B = -- B — A COB e — A' COB ^d — A" cos & d — , 

where B, A, &c., are functions of and the three other elements of its 
class. As B thus limited does not contain the three elements which are 
in the class of d, its partial differential coefficients, with respect to these 
quantities, vanish. Then the three elements accompaning in its class 
are constant, and B, as we have limited it, contains no other variables 
than and d. Thus if the differential equations determining and 6 are 
multiplied, the first by d 6, and the second bj — d 0, and the results 
added, we have an exact differential, which being integrated gives B = a. 
constant, or as it may be written, 

C= B ^ AcoBd -Jr A' cos 2 ^ + ^" cos 3 (? + 

In order to obtain the values of and 6 in terms of t, we should have 
to make another integration, but this integral suffices to show whether the 
element 0, on which depend the dimensions of the orbits, is confined be- 
tween finite limits. The value of the constant C is readily obtained by 
substituting in the right member of the equation the values of 6, and 
the three other elements of its class, which have place at a determinate 
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time, as for instance the epoch from which t is counted. Then this equa- 
tion may be regarded as the polar equation of a curve, upon which the 
values of 6 and d are always found together. Let us suppose that, d be- 
ing taken as the radius and d as the angle, the equation is represented by 
a curve having one or more 
of such branches as those 
in the figure. P\6 the pole 
from which 6 is measured, 
and P Q the line from 
which is measured. Now 
if the values of 6 and 6, 
at a determinate time, are foiind on the closed branch which envelops the 
pole P, it is plain 6 will always be comprised between certain finite limits. 
And in this case the mean motion of the argument d cannot vanish, as 
moves through the entire circumference. Here it is always possible to 
develope 6 ahd d in converging infinite series consisting of periodic terms, 
such as 

e=6, + d^ cos {e, (i + c) ] -f 02 cos 2 [»„ (^ + c) ] + 

d = d, {t+c) + d^ sin [d, (^+c) ] + d^ sin 2 [d, {t+o) ]+.... 

where 6^, 0^, 0^, &c., ^g, 6^^, 0^, &c , and o are constants. These are 
the series of which Delaunay has made such constant use in his Theory of 
the Motion of the Moon. 

But if the values of and ^ at a determinate time, are found upon the 
closed branch holding the middle place in the figure, will always be 
contained within finite limits, while d, its mean motion vanishing, will 
make oscillations forth and back betweeji definite limits. Hence although 
the mean motions are here exactly commensurable, no instability results. 
This case obtains in the three inner satellites of Jupiter, and it also has 
place in the system of the sun, earth and moon, when the last, as well as 
the earth, is regarded as a planet circulating about the sun. 

In all this we must remember that the values of and the three other 
elements of its class, both at the origin of time and ever before and after, 
must be such that they allow the development of JR in periodic series to 
be convergent; else any conclusions derived from this series are not 
legitimately established. 

From this we see that the commensurability or incommensurability of 
mean motions has no marked connection with stability. This last may be 
said to depend rather on whether the elements, such as the mean distances 
and eccentricities, which determine the dimensions of the orbits, have at 
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a given time such values as make the orbits decidedly elliptic, and permit 
to them the vibrations caused by the action of the members ot the sys- 
tem, without inteference or in other words intersections. There can hardly 
be a doubt that our solar system as composed of the sun and eight 
principal planets fulfills these conditions. 



SOLUTION OF A PROBLEM. 



BY E. W. HYDE, PROF. OF MATH. IN PA. MILITARY ACADEMY, CHESTlfiK, PA. 

Locus of a right line moving so that two fixed points of it remain con- 
stantly in two other right lines, which are perpendicular to each other but 
do not intersect. 

Let the axis of y and line -4 i? be the 
directrices. 

A line in X 3^ through two points is 

Vzz^h ^2/i-y? ui) 

/■yt fvi tyt i-fl ^ 

or, (y— 2/i)(«i-«2)=(a'-~«'i)(2/,-2/2)- 

For the line c wc have 

iBj ^ 0, y, = c sin ^, 
ajj ^ c c s 9^, 2/2 ^ 0. 

, — tan d — — . 

c cos ip X 




But 



, cos tp =^ 



r= |i_ ^'^' 



. •. sin ^ = 1/1—008^^ = 
Substituting these values in equation (a) 






or ; 



and by squaring and reducing, we obtain. 



